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The effect of fluctuations on the nuclear magnetic resonance (NMR) relaxation rate W = 
is studied in a complete phase diagram of a 2D superconductor above the upper critical field line 
Hc 2 {T) . In the region of relatively high temperatures and low magnetic fields, the relaxation rate W 
is determined by two competing effects. The first one is its decrease in result of suppression of quasi¬ 
particle density of states (DOS) due to formation of fluctuation Cooper pairs (TCP). The second 
one is a specihc, purely quantum, relaxation process of the Maki-Thompson (MT) type, which for 
low field leads to an increase of the relaxation rate. The latter describes particular fluctuation 
processes involving self-pairing of a single electron on self-intersecting trajectories of a size up to 
phase-breaking length icf, which becomes possible due to an electron spin-flip scattering event at 
a nucleus. As a result, different scenarios with either growth or decrease of the NMR relaxation 
rate are possible upon approaching the normal metal - type-II superconductor transition. The 
character of fluctuations changes along the line Hc 2 {T) from the thermal long-wavelength type in 
weak magnetic fields to the clusters of rotating FCP in fields comparable to Hc2{0). We hnd that 
below the well-dehned temperature Tq ^ O.GTco, the MT process becomes ineffective even in absence 
of intrinsic pair-breaking. The small scale of FCP rotations (^xy) in so high helds impedes formation 
of long (< i(f)) self-intersecting trajectories, causing the corresponding relaxation mechanism to lose 
its efficiency. This reduces the effect of superconducting fluctuations in the domain of high fields 
and low temperatures to just the suppression of quasi-particle DOS, analogously to the Abrikosov 
vortex phase below the Hc 2 {T) line. 

PACS numbers: 74.40.-n,74.25.nj 


I. INTRODUCTION 

Nuclear magnetic resonant (NMR) spin-lattice relax¬ 
ation is a result of nuclei-interactions with low frequency 
excitations available in the investigated system^. This 
fact makes NMR a powerful tool for studying low-energy 
excitation dynamics in novel materials^. 

In the Abrikosov phase of the type-II superconductors, 
for magnetic fields well above the critical field Hd but 
still below i^c 2 , magnetic flux lines are separated by su¬ 
perconducting regions at distances of the order of the 
coherence length ^xy The low-energy excitations driving 
spin-lattice relaxation are the weighted average of the 
intra-vortex excitations and of the contribution from the 
inter-vortex regions, possibly connected by a spin diffu¬ 
sion process^. In the vortex liquid phase, flux line dif¬ 
fusion provides an additional possible relaxation mecha¬ 
nism (see Ref. [3] and references therein). 

In a recent work^ the authors pointed out that a dy¬ 
namic state with clusters of coherently rotating FGP 
is formed above the Hd{T) line at low temperatures. 
It is therefore of special interest to study the effect of 
this fluctuation analogue of the vortex state on the mag¬ 
netic field dependence of the relaxation rate near (T). 
Some preliminary experimental studies were performed^ 
by measuring the NMR relaxation rates in a single 
crystal of superconducting YNi 2 B 2 (Tco = 15.3K in zero 
field). The authors discussed an anomalous peak in the 
NMR relaxation rate magnetic field dependence W{H) 


at temperatures 2K and 4iF in fields close to Hc 2 {T), 
which they tentatively attributed to quantum fluctua¬ 
tions of magnetic flux lines. 

Superconducting fluctuations affect the NMR spin- 
lattice relaxation rate of superconductors in a wide range 
of magnetic fields and temperatures above the upper crit¬ 
ical field line Hc 2 {T) [6-16]. First of all, they suppress 
the density of quasi-particle excitations^^’^^, which en¬ 
ters quadratically into the NMR relaxation rate, and, as 
a consequence, they reduce W. Nevertheless, this is not 
the only way for fluctuations to influence nuclear relax¬ 
ation. There exists another, purely quantum, relaxation 
process of the Maki-Thompson (MT) type which con¬ 
sists of the fluctuation self-pairing of a single electron on 
a self-intersecting trajectory due to an electron spin-flip 
scattering event at a particular nucleus^’^’^^ (see Fig. 1). 
The latter process opens a new channel of NMR signal 
relaxation leading to the increase of W. 

Below, intending to investigate first of all the region of 
the Abrikosov lattice formation from the side of normal 
metal, we study the effect of superconducting fluctua¬ 
tions both of the thermal and the quantum nature on the 
NMR relaxation mechanisms. We concentrate mainly on 
the most interesting case of a two-dimensional 5—wave 
superconductor restricting our consideration by the rep¬ 
resentative dirty limit Tcqt ^ 1, where r is the electron 
scattering time. We will derive the general expression for 
the fluctuation contribution to the NMR relaxation rate 
valid for the whole phase diagram above the line Hc 2 {T). 
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FIG. 1. (Color online) The relaxation process of Maki- 
Thompson (MT) type related to the fluctuation pairing of 
electrons on self-intersecting trajectories involving their spin- 
flip scatterings on the investigated nuclei, (a) Initially elec¬ 
tron moves along the trajectory 1 (clockwise, red), due to 
several impurity scattering events it returns to the departure 
point. In result of interaction with the nuclei, its spin and 
momentum flip and electron returns along almost the same 
trajectory 2 (counter-clockwise, green). During this motion it 
interacts with “itself in past” (corresponding superconduct¬ 
ing interactions are shaded). Such process is possible due 
to “fast” motion of the electron along the trajectory and re¬ 
tarded character of electron-phonon interaction, (b) Repre¬ 
sentation of the same process as a Feynman diagram (compare 
with Fig. 2(c). 


Its analysis for temperatures close to Tco and fields much 
less than Hc 2 { 0 ) confirms the picture of the competition 
of two contributions already studied in previous theoret¬ 
ical works^“^^. 

The situation qualitatively changes when temperature 
decreases well below Tco- The nontrivial finding con¬ 
sists of the fact that below some universal temperature 
^ 2 *D 0 — the superconducting fluctuations are no 

longer able to contribute positively to the NMR relax¬ 
ation: the coherent MT scattering is suppressed by strong 
fields H > Hc 2 {T 2 -[)), and the remaining effect of the 
quasi-particle DOS depletion results in the opening of 
a fluctuation spin gap in the magnetic field dependence 
of W. In order to compare the obtained results with 
the available low-temperature experimental data^’^^, in 


the last section we extend our theory to the case of quasi- 
two-dimensional spectrum and study the evolution of the 
crossover temperature Tq (in general, T* is a function of 
the pair-breaking parameter and we denote the lowest 
value of T* for vanishing pair-breaking as Tq*) versus the 
anisotropy parameter r = 4^^ /of a layered supercon¬ 
ductor. It turns out that three-dimensialization of the 
spectrum increases the value of T* with respect to 
which completely excludes the superconducting quantum 
fluctuations as the reason of the peak in NMR relaxation 
rate observed in Ref. [5] for T < 0.25Tco. 

The paper is organized as follows. In Sec. II we intro¬ 
duce the method of calculating superconducting fluctu¬ 
ation corrections to the NMR relaxation rate. The only 
relevant contributions to the first order in perturbation 
theory, the DOS and MT processes, are calculated sepa¬ 
rately in Sections III and IV, respectively. In Section V 
we present the total correction to the normal metal Ko- 
rringa law, which is the main result of this paper. We 
derive asymptotic expressions for the total NMR correc¬ 
tion in the regimes of quantum and thermal fluctuations 
and provide a rigorous numerical analysis of the results. 
In Section VI we present the generalization of the ap¬ 
proach to quasi-2D and 3D superconducting materials 
and outline the main consequences of the above general¬ 
ization. Finally, Section VII summarizes the main results 
of the paper and explains the physical picture behind the 
competition of the DOS and MT relaxation processes at 
different temperatures and magnetic fields in the fluctua¬ 
tion regime. The crossover temperature between the two 
regimes is obtained from qualitative considerations. 


II. MODEL 

We begin with the dynamic spin susceptibility 
X^(k, uj) = x±(k, uojy —iuj + 0+), where 

X±{k,w,) = ^'^ie“-'^(f,(5+(k,?)5_(-k,0))) . (1) 

Here S± are the spin raising and lowering operators, is 
the imaginary time, is the time ordering operator, k is 
momentum, uju = 27rTu {u = 0,1, 2...) is a bosonic Mat- 
subara frequency corresponding to the external field, and 
the angle brackets denote thermal and impurity averag¬ 
ing in the usual way. In what follows we use the system 
of units where H = ks = c = 1. The NMR relaxation 
rate W is determined by the imaginary part of the static 
limit of the dynamic spin susceptibility integrated over 
all momenta: 

A f 

W = T lim — Im / (dk) x:?_(k, —iuj ), (2) 

id J 

where (dk) = dPkj( 211 )^is the spectrum dimension¬ 
ality, H is a positive constant involving the gyromagnetic 
ratio. 
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For noninteracting electrons x^\k,cc;jy) is determined 
by the correlator of two single-electron Green’s functions 
G{k,Sn) = {iSn -^(k))“\ [Cn = + (2r)“^ sign(5n), 

= 27rT(n +1/2) is the fermionic Matsubara 
frequency, ^ (k) is the quasiparticle energy measured 
from the Fermi level], i.e. by the usual loop diagram 
with the S±(k^r) operators playing the role of external 
vertices (electron interaction with the external field). 
Its trivial calculation leads to the well-known Korringa 
law: Wq = 47rATN^{0)^ with A^(0) as the one-electron 
density of states. Below we will present the fluctua¬ 
tion contribution to W in the dimensionless form by 
normalizing it to the latter result. 

The first-order fluctuation contributions to in a 
dirty superconductor above the line Hc 2 {T) can be ex¬ 
pressed by means of the standard fluctuation “dressing” 
(see Fig. 2) of the loop of two Green’s functions by the 
fluctuation propagator Lm{^k) (wavy lines in the dia¬ 
grams) and impurity vertices Xm and Cm (shaded three- 
and four-leg blocks representing the result of averaging of 
the two Green’s functions products over elastic impurity 
scatterings in the ladder approximation). Their explicit 
expressions in the representation of Landau levels and 
Matsubara frequencies read as 


L-\^k) = 


-N{0) 



1 

2 


l^/c I + ^ 



^ 2 ) 


_T ^0(-gig2)_ 

1^1 — ^ 2 ! + co’c(m + 1/2) + T, ^ 


( 4 ) 


and the four-leg Cooperon is = [27rA^(0)r]~^ A^. 
Here m is the quantum number of the FCP Landau state, 
0/c = 2'KTk {k = 0,±1,±2...) is the bosonic Matsubara 
frequency corresponding to the FGP, 0(x) is the Heavi¬ 
side theta function. An important characteristic of these 
expressions is that they are valid even far from the crit¬ 
ical temperature [for temperatures T ^ min{r“^,cj/:)}] 
and for magnetic fields as strong as iL C iLc 2 ( 0 )/(TcO^)- 
For the sake of convenience, we introduce the reduced 
temperature t = T and reduced magnetic field 

TT^ H _ H 

where 7 ^; 1.78 is the Euler constant. The propaga¬ 

tor (3) in these variables takes the form I/“^(T, 77,0/.) = 
-N{0)£m{t,h, \k\), with 


£m (L h,x) = lnt-\- 2 p 


\ X 2h (2771 T 1) 
2 1 TT^ 


- i’ 



We will also use its derivatives £m^ (t, h, x) = 
d^£m(t^h^x)^ which can be expressed through 
polygamma functions: 




2h (2m + 1) 
t TT^ 


( 6 ) 



FIG. 2. (Color online) The diagrams for spin susceptibility. 
The solid lines correspond to free electron Green’s functions, 
bold wavy lines to the fluctuation propagator, dashed trian¬ 
gles and rectangles account for electrons scatterings on impu¬ 
rities. The two diagrams (a) represent the density of states 
(DOS) correction, the diagrams (b) represent the renormal¬ 
ization of the diffusion coefficient (DCR), while the diagram 
(c) corresponds to the Maki-Thompson (MT) process. 


Let us return to Fig. 2. The two diagrams (a) represent 
the effect of fluctuations on the single-particle self-energy, 
leading to a decrease in corresponding DOS at the Fermi 
level. Gonsequently, in accordance with the Korringa 
law, one can expect them to reduce the relaxation rate 
W with respect to its normal value, opening some kind 
of fluctuation spin-gap upon approach of the transition 
line Hc 2 (T) from the normal phase. 

Diagrams (b) with the four-leg Cooperon impurity 
blocks account for the corrections to the NMR relaxation 
rate due to the electron diffusion coefficient renormaliza¬ 
tion (DCR) by superconducting fluctuations. The anal¬ 
ogous contribution turns out to be the dominant one in 
the region of quantum fluctuations in the case of fluctu¬ 
ation conductivity^^. However, in the case under consid¬ 
eration, the additional integration over the external mo¬ 
mentum with respect to the case of conductivity makes 
their contribution proportional to the square of the small 
Ginzburg-Levanyuk number^^ 

7C(3) 

327r3N(0)Teo^2 ’ 

which strongly suppresses the entire DCR contribution^. 

Finally, the diagram in Fig. 2(c) is nothing else but the 
diagrammatic representation of the MT process shown in 
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Fig. 1, where the region of attractive interaction (in grey) 
interrupted periodically by impurity scattering events 
(circles) is replaced by the fluctuation propagator (wavy 
line). This MT type diagram for in this graphic form 
appears to be identical to the one for conductivity. Nev¬ 
ertheless, the process shown in Fig. 1(c) shows us the im¬ 
portant difference in the topology of the former and latter 
that arises from the spin structure. The MT diagram in 
Fig. 2(c) is a non-planar graph with a single fermion loop. 
In contrast, the MT graph for conductivity is planar and 
has two fermion loops. The number of loops, according 
to the rules of the diagrammatic technique^^, determines 
the sign of the contribution. In the case of spin suscepti¬ 
bility, which is under consideration, the topological sign 
of the MT diagram turns out to be opposite to that one 
for conductivity. 

The presence of the S'±(k, r) operators, taking over the 
role of external vertices, changes not only the formal sign 
of the MT diagram. The fact that two fermion lines at¬ 
tached to such vertex must have the opposite spin labels 
(up and down) eliminates the Aslamazov-Larkin diagram 
from our present consideration: one simply cannot con¬ 
sistently assign a spin label to its central fermion lines 
for spin-singlet pairing^. 


III. DOS CONTRIBUTION 


Let us start with the calculation of the DOS contribu¬ 
tion determined by the two diagrams in Fig. 2(a). The 
corresponding expression for the dynamic spin suscepti¬ 
bility integrated over all momenta is 


[ {dk) x+-^{k, = A E ^ E 

■T E 9 {Sn+u) ^k-n) , (7) 


and 






i’' 



2nT 


^(an) 


iS^k 5 ^v) 


0 ( —^/c) 0 (^/c + ^v) 



where the first summation is performed over Landau lev¬ 
els and 



(dk) G (k, Sn+v) = -i7rA^(0) sgn {Sn+u) , 

( 8 ) 


J 


DOS 


^k—n) — J (^P) ^ (P? ^n) G (p, 

N 0 (~^n^/c-n) Sgn (^n) , ^ 

= 27riA^(0)^— 2 2 • (9) 


In the approximation of a dirty metal (Tr ^ I), 

27TN‘^{0)T‘^h 


J (dk) x+?®(k,w^) 




( 10 ) 


with 


•EE Lm i^k) (U/j , Ldiy) 

m Ofc 


" /o . _^0{-en^k-n) Sgn (Sn) Sgn (en+u) 

— 2_^ O I I ’ 

(l^n ^^k—n\ T Q^m) 

( 11 ) 

where we have defined = cjc (^ + 1/2). For the 
Heaviside theta function we assume 0(0) = I. Now 
one can perform the summation over fermionic fre¬ 
quencies by splitting its domain into three inter¬ 
vals: (—oo, —— I], [—z^, —I], and[0,oo). The part 

of Eq. (10) depending on the external frequency 
which determines the imaginary part of the susceptibility 
Imx:^_(k, —iuj) in Eq. (2), is 




2 0 (^ra + Hfc) 

n=0 (^^" + Hfc + CXm) 


( 12 ) 


The summation over fermionic frequencies again can be 
performed by splitting the domain of further summation 
over the bosonic frequencies into three: k G [0,oo), k G 
[—z^, —I], and k G (—oo,—z^ —I]. Summation over the 
last interval results in zero and Eq. (12) is presented as 
the sum of regular and anomalous parts: 

+ (fi,,c..), (13) 


where 


+ 


U/j T (y.rn \ 
47rT J 




f f T \ 
V 2 ^ 47rT ) 


(14) 


^ ^ 47rT y ^ J 


(15) 


The regular part (14) is an analytic function of the exter¬ 
nal frequency uju and can be easily continued to the up¬ 


per half-plane of the complex frequencies by substitution 
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ujiy —iuj. As a result, by putting together Eqs. (2),(10), 
(11), and (14), one finds [the identity (6) was used to fi¬ 


nalize 


.(reg) 




and 


fm (^9 


Sm{l) 


(19) 


iyDOS(reg)(^^^) 


A/V(q A£'J^ {t, h, \k\) 

£m{t,h,\k\) 


(16) 


with M = {tTcor) ^ as the cut-off parameter. 



FIG. 3. (Color online) Closed integration contour C in the 
plane of complex frequencies. 


Note that the summation limit in Eq. (18) can be ex¬ 
tended from {u — 1) to z/, since fm{l = = 0. 

The analytic continuation of Eq. (18) to the upper half¬ 
plane of complex frequencies was performed in Ref. [23] 
(see also Eq. (8.84) in Ref. [19]). By means of the Eliash- 
berg transformation^"^, the corresponding sum can be 
presented as an integral over a counterclockwise closed 
contour C consisting of two horizontal lines, two vertical 
lines, and a semicircle in the upper complex plane around 
the pole 2: = 0 (see Eig. 3): 


Om (^zy) — 


1 

Yi 


coth {'Kz) — UJy)dz . 


( 20 ) 


c 

The integrals over the vertical line segments are zero and 
the integral over the semi-circle reduces to minus half 
of the residue of the integrand at ^ = 0. By inverting 
the direction of integration along the line Im z — v and 
shifting the integration variable as 2; — iujyl2'KT ^ 2; in 
the corresponding integral, one finds: 


Om (Wr.) = - - -h —y coth [WZ) 

X [fm{-iz,l^u)-fm{-iz+iVu/‘^T^T,iV^)]dz. (21) 


Now, let US proceed to the analysis of the anomalous 
DOS contribution to the relaxation rate determined by 
Eqs. (2), (10), (11), and (15). Here, the upper limit of 
summation over bosonic frequencies contains the external 
frequency and one should be cautious with the analytic 
continuation to the upper half-plane of complex frequen¬ 
cies. One can write 

f {dk) (k, = - E 

^ m 

where 


Eq. (21) is an analytic function of ujj^ and one can perform 
its continuation by the standard substitution uju —iuj. 
By the change of variables z uj j2 'kT ^ ^ in the second 
integral and with the help of the identity 

w X 1 /7N sinh(a — h) 

= -si„h(«)smh(t)- 

one finally finds 

.sinh (cc;/2T) fm{—iz,—iuj)dz 

2 J sinh {it z) sinh {it z — uj/2T) 


v-l 

Om ij^v) — ^ ^ /m(^9 
1=1 


Substitution of the explicit expression for function 
(18) fm{—iz, —iuj) from Eq. (19) into Eq. (22) results in 


j (dk) x+°®^“^(k,:^) 


N{0)h 

47r^^ 


M 


E 

m=0 


( ^r) ^m(Q) _f [^m ( ^t) ^m( \ 

Sm{ 0 ) \2T J J smh{7rz) smh{7rz — uj/2T) £m{—'^z) j 

— oo ) 

(23) 


By sending the external frequency to zero, one finds the anomalous DOS contribution to the NMR relaxation rate: 

M 




87r2^2 


m=0 


2C(0) 

£m{0) 


+ 27T 


I 


iz) lm£miiz) dz 

-oo sinh^(7r2:) \Re^£m{iz) + \m^£m{iz)\ 


(24) 
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Let us note that along the line Hc 2 {T), where the Eq. (5) 
can be simplified as 


Sm (L 


h — hc2{t) i'K^Zt 

hc2{t) 4:hc2{t) ’ 


the second term in Eq. (24) exactly cancels the first one 
and in this region only the regular part of the DOS dia¬ 
grams contributes to the NMR relaxation rate. This fact 
justifies the approximation of static fluctuations (account 
for the term with (2/^ = 0 only) assumed in the previous 
works^“^^ in their consideration performed close to Tco- 


and express Eq. (25) in the standard form 

fidk)x+I = ^ . 

^ m Ofe 

The summation over fermionic frequencies in the expres¬ 
sion 

Y' _ ^ ^/c)] ^ ['^n+z>' {^n-\-iy ^/c)] 

^ (|2£^n ^k\ T ^m) (|2^n+zy ^k\ T ^m) 

= {nk,L0,) + T^«s2+a„) 


is performed in complete analogy with the previous sec- 
IV. MAKI-THOMPSON CONTRIBUTION tion, and one finds as a result: 


The mentioned above MT contribution to the pro¬ 
cess of NMR relaxation is described by the diagram in 
Eig. 2(c). The corresponding expression for the dynamic 
spin susceptibility integrated over all momenta is 




1 1 

47rT uji, 





\^k\ 


47rT 



I ^k I T \ 
4^^T J 


(dk) x+- = y] y] Lm y] Am (e„, Qk-n) 

^ m Qfc n 

'Xmi^n+u-! ^k—n—u)^2g{^ni ^k—n)^2g{^n-\-U') ^k—n—v) ( 2 ^) 

Restricting ourselves by the assumed above case of the 
dirty superconductor, one can write an explicit expres¬ 
sion for the integral of the product of two Green’s func¬ 
tions: 

hg =J{dp)G{p, e„)G'(-p, = 27rAr(O)T0(-e„^^fe_„) 


T^eg2+a„) ^ eK 1 \y) 

AttI [LUjy + am) 




\LUy ^ am) 

/I 2LUy — \Qk\ C^m\ / / 1 

V2 ^ J 


|^/c| + 


Qlri 


47rT 


Analitic continuation of (Qki^u) is trivial, while 

that one of is performed by means of the 

Eliashberg transformation (20). Einally, one obtains 


M 




167r2^2 


m=0 

TT^ 


^ 4s:^{t,h,\k\) 

k=-oo 

oo 

dz 


+ ^(^ + 1/2) 


I, 


Im^f^(t, h, iz) 


sinn (ttz) Re^Em (L iz) + Im^f^ (t, h, iz) 


with 70 = 7r/(8TcoT0) and as the phase-breaking time. 


V. MAIN RESULT 


laxation rate, one can write the expression for valid in 
the whole phase diagram (with the restrictions discussed 
above): 


Collecting DOS and MT contributions in one expres¬ 
sion and normalizing it to the normal metal Korringa re- 


W^{t,h) 

lEo 


Gi2D f h\ 


7C(3) \ t 


TT' 


m=0 

3 


E 


.k= — oo 
oo 


X^ + f{m+l/2) 


/ 


8f^(t, h, |/c|) dz lm£'^{t,h,iz)lm£m{t,h,iz) 

- —j— 47r /- 

£m{t, h, \k\) sinh^(7r2:) Re^f^(t, h, iz) + Im^f^(t, h, iz) 

dz £m{t^h^iz) 

sinh^(7r2:) Re^f^(t, h, iz) Im^f^(t, h, iz) 


(26) 
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One can analyze it in different limiting cases. Close e = (T — Tco) ^ 1 and phase-breaking rate 7 ^ ^ 1 

to Tco and for magnetic fields not too high {h <^1) one can perform the integrations and summations in 
but arbitrary with respect to reduced temperature Eq. (26) and get 


VEfl(e,h< 1 ) 
fEo 


—3Gi2D 


{ 





168C(3) e - 7 ^ 




(27) 


In the limit of weak fields h min{e, 70 } 


of Cooper pairs is not yet strong enough: 


vpfl 


3Gi2D 


h^ 


[ ^ In ^ In 


r TT® 1 11 

. 168 C( 3 )e -70 70 e_ 

« 3Gi2D 
Wo 

_672C(3)h h_ 


70 + ^ 


168C(3) 


- 1 


(28) 


Concluding discussion of the closeness of Tco, one can 
write the explicit expression for along the line Hc 2 {T) 
in its beginning, where e + h ^ 0 : 


The first line of Eq. (28) reproduces the results of 
Refs. [6-9], while the magnetic field dependence of 
for weak fields (second line of Eq. (28)) was firstly an¬ 
alytically found in Ref. [11]. One can see that the MT 
contribution dominates when the pair-breaking is weak. 
In this case superconducting fluctuations in weak fields 
increase the NMR relaxation; increase of the field reduces 
the latter. As the phase-breaking grows, the role of the 
first term in Eq. (28) weakens and the effect of fluctua¬ 
tions can change sign: the MT trajectories shorten and 
the negative contribution of superconducting fluctuations 
due to the suppression of the quasi-particle density of 
states becomes the dominant. Since 7 ^ ^ 1 the effect of 
magnetic field on is always negative. 

In the opposite case 1 ^ h ^ max {e, 7 ^} the MT con¬ 
tribution dominates^ intrinsic pair-breaking here is 
weak while the effect of magnetic field on the motion 


f TT^ 2h , q 

W7) 1168C(3) (e + h) ( 7 ^ + h) h J ■ 

Now let us turn to the main subject of our study: the 
domain of the phase diagram above the second critical 
field at relatively low temperatures. Our general for¬ 
mula (27) allows to obtain the explicit analytical ex¬ 
pressions, for instance, along the line Hc 2 {T), where 
t hc 2 {t). Here the main contribution is due to the 
lowest Landau level of the ECP motion. Corresponding 
propagator (3) has the pole structure: 


Lo{t,h,iz) = - 


1 


1 


A^( 0 ) h ■ 




4/lc2(t) 

Performing summation over bosonic frequencies and in¬ 
tegration in Eq. (26) one finds 


w^{t « K2{t)) 
Wa 


47r2Gi2D f, 1 , 


, ( 4hc2{t)h \ 

’ [ ) 


TT^t 1 / TT^t \ 

4hc2{t)h 2 y4hc2(t)h J 


At very low temperatures t h, 
h = [H — Hc2{0)]/Hc2{0), and just above i 7 c 2 ( 0 ), 
the regime of quantum fluctuations is realized. They 
suppress the NMR relaxation due to decrease of the 
quasi-particle density of states. 


IRfl 




47r^Gi2D 1 ^ 

7C(3) h 192hl2{t)h‘^ 


(29) 


At higher temperatures, h <^t hc 2 {t), superconduct¬ 

ing fluctuations become of thermal nature, while the DOS 


suppression of the NMR relaxation remains dominant: 

K 2 (t)) _ 47r^Gi2D r 1 

'^C(3) h 16h^2(^)^ 

(30) 

The results of numerical analysis of Eq. (26) for dif¬ 
ferent pair-breaking rates are presented in Eigs. 4-5. In 
the case of a small enough pair-breaking, there is a large 
domain of the phase diagram where superconducting fluc¬ 
tuations result in the increase of the NMR relaxation rate 
(see Eig. 4). Growth of the pair-breaking suppresses MT 
contribution and when 7 ^ ^ 1 the only effect of quasi¬ 
particle DOS suppression on dominates in the whole 
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FIG. 4. (Color online) Top: The temperature and magnetic 
field dependence of the relaxation rate in case of a very 
weak pair-breaking 7 ^^ = 0.003. The thick isoline (red) repre¬ 
sents a zero relaxation rate, while the dashed isolines corre¬ 
spond to relaxation rate values of —1 and —2. The mesh-line 
t* (red) marks the critical temperature for 7 ^^ ^ 0 , while the 
light (cyan) contour line indicates the value of at hc 2 {t*) 
(—3.04) (see Fig. 6 ). Bottom: Contour plot of the region close 
to hc 2 {t*). 


phase diagram (see Fig. 5). It is interesting that even in 
the absence of the pair-breaking ( 7 ^ ^ 0 ) there exists a 
crossover temperature Tq below which the MT relaxation 
process is suppressed by strong magnetic fields and the 
fluctuation correction cannot be positive. In the case 
of a two-dimensional superconductor Tq 2 ^ ~ O.bTco* The 
temperature and field dependence of kF^(T, H) near the 
point {T 2 *Q,i^c 2 (T 2 *d)} is very singular, see the close-up 
view of its vicinity in the bottom panel of Fig. 4. 

Closer analysis of the crossover region, see Fig. 6 , re¬ 
veals that near /ic 2 (^*) the total correction curves cal¬ 
culated along the line hc 2 {t) cross at a single point 



FIG. 5. (Color online) The temperature and magnetic held 
dependence of the relaxation rate in case of a strong pair¬ 
breaking 7 (^ = 0.3. The dashed isolines correspond to relax¬ 
ation rate values of —1.5 and —2. 


at temperature ^*( 7 ^), i.e. the total correction to the 
NMR relaxation rate becomes independent of the held 
close to the point {t*,/ic 2 (^*)}- This can also be seen 
on Fig. 4, where the isoline corresponding to the value 
IT^(t*,/ic 2 (^*)) ~—3.04 is seen to be parallel to the 
t = t* line in the immediate vicinity of the superconduct¬ 
ing region. 

Below the crossover temperature t*, the total correc¬ 
tion exhibits monotonic (increasing) held-dependence for 
hxed temperature t < t*. For t hc 2 {t)^ both in the 
regime of quantum and thermal huctuations, our numer¬ 
ical analysis is in full agreement with the asymptotic ex¬ 
pressions (29) and (30), conhrming the negative sign of 
the total correction. At the same time, Fig. 6 reveals 
a non-monotonic behavior at intermediate temperatures 
t < t* when going along the hc 2 {t) line. 

Above the crossover temperature, the held dependence 
of always shows a non-monotonic behavior as a result 
of the two competing contributions, as can be seen by 
from Fig. 4. The total correction is positive (for not-too- 
strong pair-breaking 7 ^) close to the line hc 2 {t)] it then 
decreases rapidly reaching a minimum negative value at 
some intermediate distance from hc 2 {t) before increasing 
up to zero when sufhciently far from the superconducting 
region. 

Overall, our result for the total huctuation correction 
is in qualitative agreement with that of Ref. 10, 
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where the temperature range between 0.75Tco and 1.5Tco 
was analyzed to compare with experimental data. The 
authors correctly point out the strong dependence of 
on the pair-breaking parameter 7 ^. Our analysis based 
on Eq. (26) in the entire temperature range along Hc 2 {T) 
enables us to identify the temperature T*( 70 ) at which 
the DOS and MT relaxation mechanisms fully compen¬ 
sate each other, such that the fluctuation correction 
completely vanishes (in the leading order of perturba¬ 
tion theory). The dependence of this temperature on the 
pair-breaking parameter 7 ^ is presented in Fig. 6 c). The 
asymptotic crossover temperature Tq is then defined as 
T*(0), i.e. the temperature below which the negative 
DOS contribution always dominates, regardless of the 
values of 7 ^ and h. The physical picture behind this 
observation is explained in more detail in Section VII. 


VI. QUASI-TWO-DIMENSIONAL 
SUPERCONDUCTOR 

The effect of three-dimensionality of the spectrum can 
be easily accounted for by the direct generalization of the 
Eq. (26). Indeed, the properties of a quasi-two dimen¬ 
sional superconductor can be described well in the frame¬ 
work of the phenomenological Lawrence-Doniach (LD) 
model^^, which provides a Ginzburg-Landau functional 
for a layered superconductor. In the case under consid¬ 
eration, when the magnetic field is applied perpendicular 


to the layers, it takes the form: 


I 


aTcoe|^'y + ^ 




Here T/ is the order parameter of the l-th superconduct¬ 
ing layer and the phenomenological constant J is propor¬ 
tional to the energy of the Josephson coupling between 
adjacent planes. The gauge with = 0 is chosen. 

In the immediate vicinity of Tc, the LD functional 
is reduced to the GL one with the effective mass 
M = along ^-direction, where s is the inter¬ 

layer spacing. One can relate the value of J to the coher¬ 
ence length along the 2 :-direction, or, what is more con¬ 
venient, with the degree of three-dimensionality of the 
system r = 4 ^^/ 5 ^: J = oiT^^r j2 (see Ref. [19] pp. 26 
and 220). Gorresponding generalization can be done 
also in the microscopic approach: it is enough to en¬ 
rich the propagator and Cooperons (or directly the fi¬ 
nal Eq. (26)) by the account of the transversal motion: 
cjc (m + ^) ujc {rn ^ ^ ^{1 — cos Qzs) and perform 

the additional integration over the transversal momen¬ 
tum. In order to avoid the cumbersome expressions, let 
us show explicitly how it works only for the regular DOS 
contribution (16): 


VEDOS(LD,reg) 


AN{0) 

8^ V/. 


A ^ 

Ott ^ ^ 



l+|/c| I 2h (2m+l) , 4r ^*2 qzS 

2 ' t 7r2 fH 2 


Int + 7 /; 

1 + A: I 2h (2m+l) , 4r ^*2 QzS 

2 ' t 7r2 7r2t 2 



(31) 


One can see that the averaging in Eq. (31) over the 
transversal modes effectively reduces it to the same 
Eq. (16) with addition of some positive constant Xr/t 
in the arguments of all polygamma functions. Hence, in 
order to satisfy the same conditions for T* (r = 0) but at 
some temperature T* (r > 0), one should have 


hc2 (^ 2 d ) 


^2D 


hc2 (t*) 
t* 



VII. DISCUSSION 

Here we discuss the physical aspect of the results ob¬ 
tained and consequences for the general understanding 
of the fluctuation picture. As already explained above, 
it is the MT process of the self-electron pairing at the 
self-intersecting trajectories that is responsible for the 
growth of W{T,H). Its contribution to the NMR relax¬ 
ation rate is proportional to the superconducting inter¬ 
action strength g (T, H) and to the total probability for 
the formation of such trajectories (see Ref. [19]): 


This means that hc 2 (^ 20 ) /^ 2 D ^ ^c 2 (^r) A* 5 taking 

into account the monotonous increase of the second criti¬ 
cal field with the decrease of temperature, one makes sure 
that t* > ^ 20 , i.e. the temperature T* should grow with 
the increase of r. This qualitative speculation is confirmed 
by the numerical study of correspondingly generalized 
Eq. (26) (see Fig. 7), where the asymptotic crossover 
temperature is increased to tg 0.75. 


r 


inin{£^,LH}/vF 
^FCPiT,H )/vf 


X pv^dt 
Dt 


Here Lp = \/c/2eH is the FGP magnetic length, while 
^fcp(^ 5 ^) is its effective size. Glose to Tco and in weak 
fields, where the long wave-length Ginzburg-Landau fluc¬ 
tuation picture takes place = ^gl = ^xy/v^. 

In the opposite case of quantum fluctuations at zero 
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FIG. 6 . (Color online) Total correction (a) in the 2 D case for 7 ^^ = 0.003 calculated parallel to, but at different distances from 
the line hc 2 {t) in the fluctuation regime. Different colors correspond to different separations dehned by the scaling parameter 
a from hc 2 {t) given in the legend, see the (b) panel, (c) shows a linear dependence for all experimentally relevant 

values of 7 (^. For hxed temperatures below the relaxation rate is monotonically increasing with held, while for larger 
temperatures the relaxation rate is non-monotonic and grows strongly when approaching hc 2 from above. The asymptotic value 

is ti = = 0 ) = 0 . 59 . 



t 

FIG. 7. (Color online) Quasi- 2 D result for different r values 
very close to hc 2 {t) with a = 1.01 (see Fig. 6 ) for 7 ^^ = 0.003. 


temperature and in the vicinity of F 7 c 2 ( 0 ), the size 
of FCP clusters, including many pairs, is of the or¬ 
der of ^xy/V^, but the length corresponding to one of 
them is Cfcp (05 ^ Cxy (see Refs. [4] and [ 21 ]). At 
some intermediate region of temperatures along the line 
Hc 2 {T) to low temperatures the increasing magnetic held 
“breaks” GL waves and the vortex description becomes 
more adequate. The obtained crossover temperature 
T* allows us to dehne where it happens: namely for 
^FCP [T*,H{T*)] '^Lh = ^Jcl2eH{T*), where the MT 
mechanism of NMR relaxation becomes irrelevant when 
going to lower temperatures (see Figs. 4 and 6 ). The 
crossover temperature depends on the pair-breaking pa¬ 
rameter 70 and becomes minimal in the limit 7 ^ ^ 0 
with a value ~ 0.6 which can be clearly seen in Figs. 4 
and 5. 

Now let us return to discussion of the experiments of 
Ref. [5 and 20] , which partially motivated this work. The 


authors observed a well pronounced peak of W (T, H) 
versus magnetic held in the low temperature part (t = 
0 . 12 ; 0.25) when crossing the line Hc 2 {T) and attributed 
it to the possible manifestation of the quantum fluc¬ 
tuations. Unfortunately, our analysis of all fluctuation 
contributions dehnitely excludes this hypothesis: below 
to ~ 0.6 fluctuations can only open the spin gap in the 
NMR relaxation rate, but they cannot lead to its growth. 
The observed decay above Hc 2 {T) is therefore not related 
to quantum fluctuations. 
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Appendix A: Numerical evaluation of the NMR 
relaxation rate 

In order to utilize the complete expression for the NMR 
relaxation rate to analyze experimental data we need 
an efficient and accurate method to evaluate Eq. (26) nu¬ 
merically. Here we describe the method used throughout 
this work. The integral contributions (z-integrations) can 
be straight-forwardly evaluated using a suitable quadra¬ 
ture scheme. Here we use the Gauss-Legendre 5-point 
method, which also allows integration of integrable poles 
or principle values. Due to the presence of the sinh“^( 7 r 2 :) 
term in the integrand we can restrict the support to 
^ G [—5, 5]. Outside this interval the integrand is smaller 
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than the numerical accuracy (of double precision floating 
point numbers). This sum over Landau-levels is calcu¬ 
lated up to M = {tTcor)~^ explicitly. 

In contrast, the summation over k in the MT contri¬ 
bution to Eq. (26) is more involved and only slowly con¬ 
verging. For the numerical summation of the k-sum we 
separate the k = 0 term and sum from k = 1 to /Cmax 
(twice, due to symmetry) which is determined by the ar¬ 
guments of the functions being equal to = 1000. 
For k > kmax we transform the sum into an integral and 
use only the asymptotic expressions for the polygamma 
functions as the difference to the exact expression is again 
below the floating point accuracy. Then the integration 
variable is inverted and we have a finite integral for the 
remaining part of the sum. 

Therefore we concentrate on 

C^MT ^ 'Sp ft: 1^1) 


and write 


qMT 

Dm 


^max 1 


nax -L nC 

(2 —( 5 o ,/ c )+2 / 

U n kr. 


dk 


L k=0 
oMT(s) I oMT(i) 

^ m ^ m 


£m {t,h,\k\) 


QMT(i) 


= 1 Idh 

2 J Int --tp (1) +-Ip {^ + Xm) 

^max 



_ + Xm) _ 

Int - (I) - ln(2) -1- In (1 + /c + 2xm) 


with Xm = 

A convenient substitution is 


1 _ 8 {l^k)t 


1 ^ k 


z TT^ h{m-\-1/2) 7r‘^Xrj 
dz ^ ^ 

z‘^ h{m-\-1/2) 7r‘^ Xm ^ 


^max . I 2 ■ 


TT 

T 


1 kj^ 


Therefore, 


oMT(i)^![! [ dz 

8 J 


(^Y 

y TT^Xm J 


z‘^ In (t7r^x^/4) — pj (^) — ln(2) — ln( 2 :) 


h{m 1/2) 


J Am - Inz 

0 


with 


r ^ 

Ah , , 

«max = max < 2f2 — 




The sum part is calculated straightforwardly, 

which leaves the calculation of the “rest-integral” : 


with Am = In [h {m -h 1/2)] — pj (|) — ln(2). 

This integral is integrable and calculated by the Gauss- 
Legendre 5-point method (with avoids the singular point 
at ^ = 0) with only a few support points in the small 
interval 0 to Zmax using 125 support points. 

Overall this yields a highly accurate numerical value 
of the k-snms. 

In the quasi-two-dimensional case the additional finite 
g'-integral is calculated by the Gauss-Legendre 5-point 
method using 25 support points, which is sufficient to 
obtain high accuracy. 
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